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We present a theory for the van Hove singularity (VHS) in 
the tunneling density of states (TDOS) of disordered multi- 
channel quantum wires, in particular multi-wall carbon nan- 
otubes. We assume close-by gates which screen off electron- 
electron interactions. Diagrammatic perturbation theory 
within the non-crossing approximation yields analytical ex- 
pressions governing the disorder-induced broadening and shift 
of VHS's as new subbands are opened. This problem is non- 
trivial because the (lowest-order) Born approximation breaks 
down close to the VHS. Interestingly, compared to the bulk 
case, the boundary TDOS shows drastically altered VHS, even 
in the clean limit. 



Van Hove singularities (VHS's) in the thermodynami|C| 
density of states (DOS) have been predicted in 1953EJ 
and were observed in many experiments since then. The 
DOS for a (i-dimensional system with dispersion relation 
E(k) can be written as an integral over the Fermi surface, 
v(E) = (2ir)- d f dS/\V%E(k)\. The quantity in the de- 
nominator is basically the group velocity. Due to symme- 
tries in a crystal, the group velocity may vanish at certain 
momenta, resulting in a divergent integrand. This di- 
vergence is integrable in three dimensions, and typically 
leads to a finite DOS. In lower dimensions, however, very 
pronounced VHS appear. In this paper, we focus on the 
one-dimensional (ID) limit, where the VHS diverges like 
v{E) ~ l/yE — E n in a clean system when the energy 
E approaches the threshold E n from above. Therefore 
VHS's appear as sharp features in the energy-dependent 
DOS reflecting the onset of new active subbands. Similar 
VHS's exist for the tunneling density of states (TDOS) 
measured at some location x along the system, 
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v(E,x) = ^ / dte iEt l h {i>(x,mUxM, 

where ip(x, t) is the electron operator. The TDOS is eas- 
ily accessible experimentally by measuring the conduc- 
tance through a weak link or a tunnel junction. Typi- 
cally, it is obtained by means of scanning tunneling spec- 
troscopy (STS).Ert2 In general, one has to carefully dis- 
tinguish between the DOS and the TDOS. While the 
DOS is a thermodynamic property of the system, the 
TDOS is a local property and therefore depends on the 
position. In the "bulk" limit, the DOS and the (possi- 
bly coarse-grained) TDOS are expected to be identical, 
but near a boundary they can strongly differ. In the 
present paper we address the question: How is the VHS 



modified by disorder? This question is of importance 
for the interpretation of experiments on multi-wall nan- 
otubes (MWNTs). A MWNT consists of a few concentri- 
cally arranged graphite shells, with a typical outermost- 
shell radius R w 5 to 20 nm and lengths of up to 1 mm. 
As demonstrated by magnetoresistance measurements, □ 
electronic transport occurs only through the outermost 
shell and exhibits the typical fingerprints of diffusive be- 
havior while single-jarall nanotubes (SWNTs) are ballistic 
ID quantum wiresfl The observation of VHS's in recent 
STS experiments on SWNTs on a metallic substrate -rep- 
resents a direct proof for a ID band structure.BOfi In 
particular, the predicted 1 / \JE — E n behavior of ballis- 
tic ID wires has been observed. 

Since in MWNTs the incommensurate inner-shell ionic 
potential acts on electrons in the outermost shell, the re- 
sulting quasi-periodic potential effectively acts as a ran- 
dom potentialE We then want to understand how VHS's 
develop with increasing amounts of disorder. By inten- 
tionally damaging MWNTs, e.g. by fast ion bombarding, 
any level of disorder may be realized experimentally to 
test our predictions below. For typical MWNTs with 
mean free path I « 2irR, the characteristic subband fea- 
tures of the TDOS should still be present, albeit consid- 
erably broadened and possibly shifted. This question is 
clearly also of relevance to other quasi-lD quantum wires 
such as long chain molecules. Previous experiments on 
intrinsically, (hole-) doped MWNTs reported by Bach- 
told et alS3 have observed power-law zero bias anoma- 
lies due to electron-electron interactionsB Here we focus 
on a completely different scenario, where interactions are 
screened off by working on a metallic substrate which 
is typical for STS measurements. Hence we treat only 
the non-interacting problem. Then spin only contributes 
trivial factors of two and can be ignored. In addition, 
disorder-induced scattering between the two distinct gap- 
less Fermi (K) points of the first Brillouin zone of the 
honeycomb lattice is expected to be largely suppressed, 
and we thus consider only one Fermi point. Since in typ- 
ical STS experiments only the TDOS of the outer-most 
shell is probed, in the following an effective single-shell 
model is assumed, where inner shells only give rise to a 
disorder potential for outer-shell electrons. 

Nanotubes (NTs) can be thought of as graphite sheets 
wrapped onto a cylinder. The low-energy theory of a 
clean graphite sheet is given by the 2D Dirac Hamilto- 
nian. This graphite sheet is rolled up into a tube by 
enforcing periodic boundary conditions around the y di- 
rection, the tube pointing along the x direction. For later 
convenience, we define = Tivf/2'kR with vf = 8xl0 5 



1 



m/s. (To simplify notation, we often set Ti — vf 
the absence of disorder, we have 
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with k = (k,uj n ). The "Matsubara frequencies" uj n given 
by E n = 2im/j3 = hu> n with integer n arise due to the 
finite radius, ip is a two-component spinor, where the two 
components reflect the sublattice degree of freedom, as 
the honeycomb lattice has a basis containing two atoms, 
and the Pauli matrices a — (cr x , a y ) act in this space. The 
"mass" M is generally nonzero for NTs due to chirality 
effecte,or an applied magnetic field B parallel to the tube 
axis .til For a flux <I> (in units of the flux quantum), we 
get M = 27r<£>//3. Since chirality effects can always be 
absorbed by an adjustment of B, we take $ = B/Bq 
with Bo = h/{eirR 2 ). The mass term in Eq. (|l|) couples 
in the same way as k y = u) n , and hence we can replace 



uj n -> u n + 2tt$//3 = (n + B/B )/R 
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in Eq. (Q) to include the mass term. For integer <E>, 
the system is obviously not influenced by the magnetic 
field, as we can simply shift n to absorb <£>. Next we 
discuss our modelling of disorder. On the one hand, 
there may be disorder-induced hopping events between 
nearest-neighbor sites on the honeycomb lattice via real 
or virtual states provided by the inner-shell ionic poten- 
tial. Since hopping connects different sublattices, the re- 
sulting modulation of the hopping matrix element leads 
to a random gauge field. Gauge field disorder also cap- 
tureS|-thc effects of topological defects in the outermost 
shell.ll 2 ] However, in particular once the MWNT is in- 
trinsically damaged, the dominant disorder mechanism 
should be due to direct impurity potential scattering pro- 
cesses which are diagonal in sublattice space, 
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In what follows, we retain only disorder of the type 
(|J). The scattering potential is taken as a static Gaus- 
sian random field with zero mean and (V(x)V(x')) — 
Ay<5(x — a?), where Ay = h 2 v F A/R defines the dimen- 
sionless disorder strength A. 

In addition to the bulk case, we also address the bound- 
ary TDOS arising when one tunnels into the end of a 
MWNT. In reality, the end TDOS could of course be 
quite complex due to the formation of bound states. A 
few lattice spacings away from the end, however, we ex- 
pect that the situation can be described by a continuum 
model. Surprisingly, the presence of a boundary implies 
a drastic change in the TDOS even in the absence of 
interactions, namely a strong suppression of the TDOS 
close to the boundary. Since this can be demonstrated al- 
ready in the absence of disorder, let us briefly discuss the 
clean limit. The TDOS at position x and energy E > 
is v{E,x) = —IiaTv^ iy G(E,k,x)/'!T, where we have to 



trace over the "spin" and the momentum. We consider a 
semi-infinite [x > 0) tube, assuming a hard-wall poten- 
tial at the boundary x — 0. The Greens function is then 
G (E, k,x) = 2 sin 2 (kx) / (E - a ■ k + i0 + ), and hence 
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where is the Heaviside step function and vid — I/ttHvf 
serves as natural unit for the TDOS (without spin and K 
point degeneracy). As one can see from Eq. (|J) with (0), 
a gap is generally present. This gap varies as a function of 
B, and vanishes periodically with period Bq. Far away 
from the boundary, sin 2 (x^_E 2 — E 2 ) — » 1/2, and we 
obtain the bulk TDOS (see also Ref. JL|) 
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which equals the thermodynamic DOS. The well-known 
1 / \JE — E n VHS's in ID appear at the onset of new sub- 
bands. Since the summation also includes negative values 
of n, a magnetic field in general causes a doubling of the 
VHS's, see Eq. (||). Sufficiently close to the boundary, 
however, we obtain a completely different result. From 
Eq. (|) we find the boundary TDOS 
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For E — ► 0, this predicts v(E) ~ E 2 for B — 0, in con- 
trast to the finite boundary TDOS for a doped tube. This 
behavior can be traced back to the linear dispersion re- 
lation of Dirac fermions. More interestingly, the typical 
ID VHS of the bulk TDOS is drastically altered close to 
the boundary. Instead of a divergence, the only sign of 
the opening of new subbands is a nonanalyticity at the 
threshold energy E n , with a square-root energy depen- 
dence of the boundary TDOS above E n . This possibly 
explaines the observation of Ref. ||, where no divergen- 
cies were found in the tunneling spectrum at the end of 
a SWNT. We note that this phenomenon is quite gen- 
eral. For both Dirac and Schrodinger fermions, the ex- 
ponent in the TDOS energy dependence close to a VHS 
changes by one when going from the bulk to the bound- 
ary limit. For E — > 0, however, for Dirac fermions, the 
exponent changes by two due to the special role of the 
energy E = 0. The relevant crossover scale between bulk 
and boundary behavior of the TDOS depends on energy. 
Focussing on E close to but above a given thr eshold en- 
ergy E n , this scale is x* « hvp/^/E 2 — E 2 t . At zero 
temperature, the bulk limit is reached for x 3> x* , and 
the boundary limit for x <C x* . For finite temperatures 
T, if the thermal scale xt = TiVF/ksT is smaller than 
x*, one should replace x* by xt- 
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We study disordered MWNTs using diagrammatic per- 
turbation, theory within the non-crossing approximation 
(NCA) .til The use of NCA for this problem requires some 
care. On the one hand, for 2D random Dirac fcrmions, di- 
agrams with crossing impurity lines cause the same (log- 
arithmic) singularities as rainbow diagrams in arbitrary 
orders of perturbation theoryjlj and crossing diagrams 
must be treated on the same footing as the rainbow ones. 
On the other hand, in the single-channel ID limit, it is 
also well-known from exact calculations that crossed di- 
agrams are in general as important as rainbow diagrams. 
Since MWNTs are in between the ID and the 2D limit, 
studying the importance of crossed diagrams is there- 
fore mandatory. The situation for MWNTs is fortunately 
quite different from both limits mentioned above. First, 
there are no logarithmic singularities appearing in the 
self energy expansion as in the 2D limit. Second, evalua- 
tion of the simplest diagrams using standard methodsEil 
shows that rainbow diagrams are larger than crossed di- 
agrams by a factor A-kv[E)/v\d- For the energy range 
of interest here, we are therefore entitled to compute the 
TDOS within NCA. We focus on the bulk case here and 
briefly comment on the boundary result later. 

With G (E, k) = (E - a ■ k + i0+ ) - 1 , the self energy to 
lowest order in Ay (Born approximation) takes the form 
= -in A v vq( e )> with v (E) given by Eq. (§). 
Remarkably, the Born approximation breaks down in the 
vicinity of a VHS. Specifically, the second-order self en- 
ergy contribution obeys 
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which diverges for E approaching E n from above. We 
therefore must address the higher-order contributions to 
the self energy. Even within NCA, since the perturba- 
tion expansion is asymptotic, we have to arrange the 
order of summation in a physically meaningful way to 
avoid familiar but unphysical divergencies and incon- 
sistencies. This would also be important for analyzing 
diagrams beyond NCA, see Ref. To achieve that, 
we follow the self-consistent iterative approach proposed 
by LeeJl3 Within this approach, the self energy Ejy in- 
cluding all contributions up to Ath order (N > 1) is 
Y,jy(E) = AyTrg a GN-i(E, k), with the corresponding 

Dyson equation G]^(E,k) = G Q 1 (E,k) - T, N (E). This 
form is then used to calculate the self energy Ejv+i. In 
the limit N — > oo, this procedure converges leading to 
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which has to be solved self-consistently for E(_E). One 
can check explicitly in each order that Eq. (^) reproduces 
all non-crossing diagrams correctly. The result can then 
be used to compute the TDOS, v{E) = -Im E(E)/ttA v . 
It is not possible to simply assume an energy-independent 
mean free path for all energies. The energy-dependent 



mean free time is — Im£(i?) = h/2r(E), and therefore 
also an energy-dependent mean free path 1(E) — vft(E) 
results. 

From Eq. (0), we obtain (E R = Re £) 
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where correct units were restored for clarity. Equation 
(||) can easily be solved numerically. Once we know the 
self energy, the TDOS follows as 



v{E) = Rc ^. (E - £(£)) sgn(E - E H ) 



^{E-H(E)Y-El 



(9) 



For A > 0, the finite imaginary part of E in the denomi- 
nator of (||) causes a broadening of the VHS, whereas the 
real part causes an energy shift of the peaks. For the nu- 
merical evaluation of Eq. (||), a cutoff for the band index 
n has to be used, which is naturally given from the band- 
structure. For instance, for armchair NTs, the number of 
subbands is given by 2A = 8irR/V3aM For R = 10 nm, 
we get N w 295. The results are, however, not very 
sensitive to the precise choice of this cutoff. Equations 
(||) and (0) can be used to fit experimental data for the 
TDOS of MWNTs. Assuming that R is known, since A is 
the only fit parameter, the disorder strength can then be 
determined directly from the TDOS which should allow 
for detailed comparison to STS experiments on MWNTs. 
This would provide precious information on the disorder 
strength. 

Figure [l] shows the strong broadening of the VHS due 
to disorder. For sufficiently strong disorder, namely when 
motion around the circumference becomes diffusive, the 
peaks can even disappear completely above a certain en- 
ergy which decreases with increasing A. In the region 
where no VHS's are present, the TDOS apparently be- 
haves like a power law, v(E) ~ E a , with a(A) < 1 This 
power law scaling holds remarkably well for A < 0.05. 
Note that energies are not measured relative to Ep, and 
therefore this power law is unrelated to the findings of 
Rcf. [l(]. Also it has a different origin than the power 
law found in the 2D-case that rests upon the inclusion 
of crossed diagrams.Ej In addition, the positions of the 
VHS's are shifted to smaller energies with increasing A. 
This shift grows approximately linearly with A, and the 
relative shift, compared to the position of the VHS in the 
clean system, can be up to 20% depending on the disorder 
strength. Since the radius of a NT is often determined 
from the relative positions of the VHS, this observation 
suggests that such interpretations need to be taken with 
some caution. The disorder-induced shift has to be taken 
into account to obtain correct results. With decreasing 
radius and for increasing order n of the VHS, the rel- 
ative shift becomes systematically smaller. Next, if we 
consider a fixed A and vary a magnetic field B applied 
parallel to the tube axis, we find the situation depicted in 
Fig. ||. Again, there is a doubling of the VHS's, and the 
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shift of the positions varies periodically with B. As in the 
clean case, there is a gap even in the presence of disorder. 
But due to the A-dependent shift of the VHS towards 
smaller energies, the gap gets partially filled with states 
and is therefore smaller than in the clean case. How- 
ever, magnetic field effects are only weakly affected by 
disorder. Finally, we briefly turn to the boundary limit. 
Since it is easy to generalize Eqs. (||) and @ to obtain 
the boundary TDOS, we only mention the modifications 
compared to the clean case. The disordered boundary 
TDOS is suppressed compared to the clean case, where 
the suppression increases with increasing A. The posi- 
tions of the van Hove nonanalyticities at the opening of 
new subbands are shifted to higher energies, although the 
shift in the bulk case was to smaller energies. The form 
of the nonanalyticity is not significantly changed by dis- 
order, but approximately retains the square-root energy 
dependence. 

To conclude, we have calculated the TDOS of disor- 
dered multichannel quantum wires, with special empha- 
sis on MWNTs. In the present theory, electron-electron 
interactions are supposed to be screened off by a metallic 
substrate or a close-by gate. Focussing on static disorder, 
within NCA a self-consistent non-pertubative summation 
of all diagrams for the self energy yields an analytical 
result for the disorder-broadened VHS's. For given ra- 
dius of the MWNT, this result involves only one param- 
eter (the disorder strength A), which should allow for 
a detailed comparison to STS experiments on MWNTs, 
where the disorder strength can be tuned at will. 
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FIG. 1. Bulk TDOS of a disordered tube with R = 10 nm 
at B = for different values of A. The broadening of the 
VHS with increasing A is clearly visible, as well as their shift 
towards smaller energies. 
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FIG. 2. Bulk TDOS for A = 0.01, R = 10 nm and different 
values of B. The curves corresponding to different B are 
shifted vertically by the same amount for better visibility. 
The lowest curve corresponds to B = 0, then B increases in 
steps of AB — 2.2 T. Notice the opening and closing of a gap 
with period B = 13.2 T. 
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